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Abstract. The aim of this paper is to determine the initial rotational ve- 
locities required on the ZAMS for single stars to reach the critical velocity, 
sometimes called the break-up velocity, during the Main-Sequence (MS) phase. 
Some useful relations between f2/f2 cr ;t, f/w C rit (w is the velocity at the equator), 
the moments of inertia, the angular momenta, the kinetic energy in the rotation 
and various other basic physical quantities are obtained. 



1. Stellar evolution near the critical limit 

The n ew results from VLTI on the shape of Achernar ((Domiciano de Souza et all 
2005) show that detailed parameters of rotating stars become accessible to in- 
terferometric observations. In order to make valuable comparisons with models 
of rotating stars, we provide here a few basic data on rotating stellar models. 

We also search the initial conditions required for a star to reach the crit- 
ical limit during its Main-Sequence phase. When the surface velocity of the 
star reaches the critical velocity {i.e. the velocity such that the centrifugal 
acceleration ex actly balances gravity ), one expects that equatorial ejections of 
matters ensue ( Townsend et al.1l20 04b The physics involved in the ejection pro- 



cess, determining the quantity of mass ejected, the timescales between outburst 
episodes, the conditions required for a star to present such outbursts are impor- 
tant questions not only for understanding Be, B[e] or Luminous Blue Variable 
stars, but probably also for having a better knowledge on how the first stellar 
generations in the Universe behaved ((Ekstrom et al.ll2005h . Indeed, at lower 
metallicity, the stellar winds are weaker and much less angular momentum is 
removed at the surface. This favors an evolution to the critical limit and may 
have important consequences for the evolution of the first stellar generations. 
Also there are some indications that the distrib ution of the initial velocities 
contains more rapid rotators at lower metallicity (M aeder et alJll9 99). Finally, 
realistic simulations of the formation of the first stars in the Universe show that 
the problem of the dissipation of the angular momentum is more severe at very 
low Z than at the solar Z. Thus these stars might begi n their evolution with a 
higher amount of angular momentum ( Abel et al.ll2002l ). 



Rotation affects all the output s of the stellar models. The reader will find 
reviews on the effects of rotation in iMaeder fc Mevnetl (l2000aT l . 15 eger fc Langer 



: 2000b I'l aloul (1200 ll ) and lMevnet fc Maederl (l2005bl b In this work, we compute 



112 different stellar models with initial masses equal to 3, 9, 20 and 60 M , for 
metallicities Z equal to 0, 0.00001, 0.002 and 0.020, and for values of the ratio 
of the angular velocity to the critical angular velocity fi C rit equal to 0.1, 0.3, 
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Figure 1. Left: Variations of the momentum of inertia on the ZAMS as a 
function of the initial masses for various metallicities Z (squares: Z = 0.020; 
triangles: Z = 0.002; pentagons: Z = 0.00001 and circles: Z = 0). Right: 
Variations of the gravitational energy on the ZAMS as a function of Mini for 
various Z. 

0.5, 0.7, 0.8, 0.9, 0.99. For the 3, 9 and 20 M models, the computation was 
performed until the end of the Main-Sequence phase or until the star reaches the 
critical velocity. For the 60 M stars, we have obtained all the models on the 
ZAMS, but only a subset of them were computed further on the Main-Sequence. 
We consider that a star arrives on the ZAMS, when a fraction of 0.003 in mass 
of hydrogen has been burned at the center. On the ZAMS, the star is supposed 
to have a solid body rotation. During the Main-Sequence phase the variation of 
Q inside the star is computed self-consistently taking into account the various 
transport mechanisms, i.e. convection, shear diffusion, meridional circulation 
and horizontal turbulence. The removal of angular momentum at the surface by 
the stellar winds and the changes of f2 resulting from movements of contraction 
or expansion are also accounted for. A detailed description of the physical 
ingredients used in this grid of models will be given elsewhere (Ekstrom et al., 
in preparation), let us just mention here that the prescriptions for the opacity 
ta bles, the initial chemical co mpositions and the treatment of convection are as 
in lMevnet fc Maederl (|2005aT l. 

Based on the set of numerical results described above, we first present use- 
ful relations between basic quantities, and then explore the initial conditions 
required for a star to reach the critical velocity during it MS evolution. 

2. Moments of inertia and energy contents on the ZAMS 

Figure^shows for different initial masses and metallicities the moment of inertia 
and the gravitational energy of the stellar models on the ZAMS. Only the results 
obtained for a value of Q/Q CT i t equal to 0.7 are presented, since these quantities 
present a very weak dependence on Q/0, CT n. Numerically, these quantities were 
computed by summing the contribution of the different shells composing the 
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Figure 2. Left: Variations of the polar radius on the ZAMS as a function 
of Mi n i for various metallicities. Right: Fraction of the total mass (m r is the 
mass inside the radius r) as a function of the fraction of the total radius for 
solar metallicity models with J7/f2 cr it=0.1. 
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Figure 3. Left: Variation of the critical equatorial velocity on the ZAMS 
as a function of the initial mass for various metallicities. Right: Variation of 
the critical rotation period on the ZAMS as a function of the initial mass for 
various metallicities. 
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Figure 4. Variations of the angular momentum on the ZAMS as a function 
of f2/fi cr it for various metallicities, plotted according to the initial mass of the 
models. Symbols as in Fig. ^ The angular momentum is given in units of 
10 53 g cm 2 s _1 . Same symbols as in Fig. ^ 
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Figure 5. Ratio of the kinetic energy in the rotation to the gravitational 
energy as a function of Q/tt CT it for various metallicities, plotted according to 
the initial mass of the models. Same symbols as in Fig.^ 



star (between 250 and 400 shells). Since the deformation due to rotation is 
very weak in most of the star (even near the critical velocity) and since the 
outer layers contain little mass, we neglect, in all the following estimates, the 
deformation of the shells induced by rotation. To each isobaric surface, limiting 
a volume V (which is not a sphere), we associate an average radius r given by 
r=(F/(4/3» 1 /3. 

Looking at the results concerning the moments of inertia, we note that, at 
a given metallicity, the moments of inertia increase with the initial mass. The 
increase amounts to about 3 orders of magnitudes between 3 and 60 Mq. For 
the non-zero metallicities, we have linear relations between log I and log M. 
In the case of the Pop III stellar models, the linear relation breaks down and 
the increase of the moment of inertia with the mass is less steep than at higher 
metallicities. 
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What is the cause of this difference ? Passing from 3 to 60 M increases the 
radius by a factor 5-6, whatever the metallicity between Z = 0.00001 and 0.020 
(see the left part of Fig. On the other hand, for zero metallicity stars, the 
radius passes from about 1 R for a 3 M stellar model to 3.2 R for a 60 M , 
thus the increase amounts to only a factor 3. This comes from the fact that, at 
the beginning of the core H-burning phase in Pop III massive stars, only the pp 
chains are active. The pp chains are not efficient enough to compensate for the 
energy lost at the surface and the star must extract energy from its gravitational 
reservoir. The contraction is more severe in the higher mass range where the 
luminosities are higher. This explains the shallower increase of the radius as a 
function of the mass in Pop III stars and therefore the less steep increase of the 
moment of inertia with the mass. 

Let us add here that the contraction occurring in Pop III stars at the begin- 
ning of the core H-burning phase leads to an increase of the central temperature 
until the point of helium ignition is reached. Then, due to He-burning, small 
amounts of carbon and oxygen are synthesized. When the abundances of carbon 
and oxygen in the core reach a level of 10~ 10 in mass fraction, the H-burning can 
occur through the CNO cycle and this cycle becomes, as at higher metallicity, 
the dominant energy source for the rest of the Main-Sequence phase. 

Looking at the left part of Fig. ^ we note that, for a given initial mass, the 
variation as a function of the metallicity is much smaller than the variation as 
a function of the initial mass at a given metallicity. There is at most an order 
of magnitude of decrease when the metallicity passes from Z=0.020 to 0. The 
decrease is more pronounced in the range of the high mass stars than in the 
range of the intermediate mass stars. The reason for that is again the difference 
in the way the beginning of the H-burning occurs in massive Pop III stars. 

The gravitational energy or binding energy is shown in the right part of 
Fig. ^ Passing from 3 to 60 M , the binding energy increases by about two 
orders of magnitude. For a given initial mass, the dependence on the metallicity 
is much more modest. For a 20 M stellar model at Z = 0.020 the binding 
energy is about 4-10 50 ergs, its Z = counterparts has a binding energy equal 
to slightly more than 10 51 ergs, i.e. about a factor 2.5 higher. Let us recall for 
comparison that the binding energy of a neutron star is of the order of 10 53 ergs, 
about two orders of magnitude greater. 

The left part of Fig. |2]shows how the polar radii of the different models vary 
as a function of the initial mass and metallicity. The right part of this figure 
illustrates the fact that in the 3 M stellar model, a given fraction of the total 
mass is enclosed in a smaller fraction of the total radius than in the 60 M , 
indicating that lower mass stars have a more concentrated mass distribution 
than higher mass stars. 

3. Critical velocity, angular momentum, kinetic energy of rotation 
and surface velocity on the ZAMS 

The variation of the critical velocity as a function of the initial mass and metal- 
licity is shown in Fig. At a given metallicity, the critical velocity is higher 
for higher initial mass stars. The increase is not far from being linear. The 
slope remains nearly the same for metallicities higher than 0.00001, while it is 




Figure 6. Variations of the equatorial velocity at the surface on the ZAMS 
as a function of f2/f2 cr it for various metallicities, plotted according to the 
initial mass of the models. Same symbols as in Fig. ^ 
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Figure 7. Left: Variation of u/u C rit as a function of u> = f2/f2 cr ;t. Right: 
Variation of the ratio Re(oj) / Rp(co) (equatorial over polar radius) as a function 
of oj = n/Jlcrit- These relations have been obtained in the frame of the Roche 
model. They do not depend on the mass, the metallicity or the evolutionary 
stage considered. 

much steeper for the pop III stellar models. For a given initial mass, the critical 
velocity is higher at lower Z. For the 60 M model, the critical velocity passes 
from slightly more than 800 km s _1 at Z=0.02 to nearly 1600 km s _1 at Z=0. 

The critical periods are shown on the right part of Fig. |3J They are related 
to f2 cr it by the simple relation P cr jt = 2-7r/f2 cr it. The critical periods are com- 
prised between 2.5 and 25 days for the whole domain of masses and metallicities 
explored in the present work. 

The variation with f2 /fi c rit of the total angular momentum content is shown 
in Fig. @J Note that the highest point on the right of the figures gives an esti- 
mate of what could be called the critical angular momentum content. Since the 
moment of inertia of the star does not depend much on the rotation velocity, 
and since on the ZAMS, we supposed solid body rotation, one obtains linear 
relations between C and S7/f2 cr i t . We note that the metallicity dependence re- 
mains modest at least in the range of metallicities between 0.00001 and 0.020. 
For the 9 M models and above, starting from the Z=0.020 models, there is 
first a slight increase of the angular momentum content when lower metallicities 
are considered, the initial mass and the value of fi/J7 cr it being kept constant 
(compare the curve with the black squares for the Z=0.020 models with the 
curve with the triangles corresponding to 0.002 models). Then for still lower 
metallicities, the angular momentum decreases with the metallicity (the curve 
for the Z=0. 00001 models overlap the Z=0.02 models) reaching its lower values 
for the Pop III stellar models. 

In the mass range between 9 and 20 M , two stars on the ZAMS, having 
the same angular momentum content, would have very similar value of Q/0, CT [ t 
whatever their metallicity between 0.00001 and 0.02. Only if the metallicity 
is zero, would the value of S7/f2 cr it corresponding to that angular momentum 
content be much higher. As a numerical example, a 20 M on the ZAMS with £= 
0.3T0 53 cm 2 g sec -1 has a Q/Q CT i t = 0.5 when Z is comprised between 0.00001 
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Figure 8. Left: Variation of the polar radius R p (u>) as a function of u> = 
f2/f2 cr it, normalized to the value of R p at w = 0.10, for various initial masses 
at Z = 0.020. Right: Variation of the total luminosity L(uS) as a function of 
ui = n/ficrjt), normalized to the value of L at u) — 0.10, for the 60 Mq models 
at various metallicities (squares: Z = 0.020; triangles: Z = 0.002; pentagons: 
Z = 1Q- 5 and circles: Z = 0). 



Table 1. Velocity in km s 1 on the ZAMS corresponding to f2/f2 cr ;t = 0.7 
M Z=0 Z=0.00001 Z=0.002 Z=0.020 



3 312 301 262 224 

9 473 383 325 287 

20 644 447 380 344 

60 795 566 480 428 



and 0.020. This angular momentum content corresponds to f2/f2 cr it = 0.8 when 
Z=0. 

The rotational kinetic energy expressed as a fraction of the binding energy 
is shown in Fig. |21 Rotational energy amounts to at most a percent of the 
binding energy. This is consistent with the well known fact, that the effects of 
the centrifugal acceleration remains quite modest on the hydrostatic structure 
of the stellar interior, even when the surface rotates with a velocity near the 
critical one. 

The relations between the equatorial velocity and f2/$7 cr it are plotted in 
Fig. El We give in Tabled the values of v obtained for $7 / S7 C rit = 0.7 for various 
masses and metallicities. The value £l/Q cr [ t = 0.7 corresponds to an initial 
rotation of 290 km s _1 for a 9 Mq stellar model. The time averaged velocity of 
this model during the Main-Sequence phase is equal to ~240 km s _1 . 

We see that for a given value of r2/f2 cr it and of the initial mass, the surface 
equatorial velocity is higher at lower Z. At a given metallicity, a value of $l/£l CT [ t 
corresponds to a value of the equatorial velocity which is higher for higher initial 
masses. This means that if the stars begin their evolution on the ZAMS with 
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approximately the same value of f2/f2 cr i t , they have a higher surface velocity at 
lower metallicities and for higher initial masses. For instance, a 3 M stellar 
model with an initial value of fi/O cr i t equal to 0.7 has an initial velocity equal 
to 220 km s _1 (see Table ^), while a 60 M Q model, starting with the same value 
of f2/$7 cr j t , has an initial velocity of 800 km s _1 . 



4. The Roche model and surface deformation on the ZAMS 

The relation between f/f cr it and Q/Q CT i t obtained in the frame of the Roche 
model is shown in Fig. [7| (left part). This relation is the same whatever the 
initial mass or metallicity. One sees that except at the two extremes, the values 
of v j v cr it is smaller that that of Vt/ Sl C rit • The variation with 0, /fl cr it of the ratio 
of the equatorial radius to the polar radius is also a property of the Roche model. 
It is shown in the right part of Fig. [7| We note that for f2/fj cr j t below 0.7, the 
equatorial radius is longer than the polar one by less than 10%. 

The variation of the polar radius when the rotation velocity varies is illus- 
trated in the left part of Fig. |HJ For the 3, 9 and 20 M at solar metallicity, 
the polar radius decreases when the rotation increases. The decrease is at most 
1.5%. For the 60 Mq model the behaviour is reversed, the polar radius slightly 
increases when rotation increases. The maximum increase amounts to 4.5%. 
This difference of behavior for the 60 Mq star is likely linked to the fact that in 
massive stars radiation pressure plays a more important role. 

In the right part of Fig. |S1 the variation of the luminosity as a function 
of the rotation is shown for 60 M Q stellar models on the ZAMS at various 
metallicities. One sees that the luminosity decreases by at most 5 to 10% when 
rotation increases up to the critical limit. This is due to the fact that a rotating 
star is, on the average, more extended. As a consequence, the gradients of the 
temperature and thus the fluxes are smaller. 



5. Initial conditions for reaching the critical limit 

The conditions for reachi ng the break-up limit have recently been reviewed by 
iMevnet fc Maederl ^006). We concentrate here on the case of 20 M stellar 



models. The cases of the other initial masses will be discussed in Ekstrbm et al. 
(in preparation). 

Figure shows the evolution of the surface velocity expressed as a fraction 
of the critical velocity during the Main-Sequence phase. The computations were 
stopped as soon as the star reaches the critical limit. Quite generally the curves 
can be decomposed into three parts, whose importance depends on the initial 
rotation and metallicity: 

1. At the very beginning, there is a short adjustment period, which lasts 
for a few percents of the Main-Sequence lifetime, during which meridional 
circulation transports angular momentum from the outer parts of the star 
to the in ner ones, u ntil an asymptotic profile of the angular velocity is 
reached (jZahnl H992). As a consequence the surface velocity decreases. 
For instance, a value of v/v cr i t on the ZAMS equal to 0.35 passes to a 
value of 0.30 or slightly lower after a few 10 5 years. For lower initial values 
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Figure 9. Evolution as a function of time of v/v C rit at the surface of 20 M Q 
stellar models for different initial metallicities and velocities. From bottom to 
top, the models with initial values of f2/r2 C rit equal to 0.1, 0.3, 0.5, 0.7, 0.8, 
0.9 and 0.99 are shown. 



of v/v cr it the decrease is much weaker, while for higher values, it is higher. 
This illustrates the dependence of the meridional velocity on f2. 



As explained for instance in lMaeder Mevnetl ( 2000al ). after the first ad- 
justment period, a large outer cell of meridional circulation sets in, trans- 
porting angular momentum from the inner parts of the star to the outer 
ones. If this transport is rapid enough and the stellar winds not too in- 
tense, then, v/v cr - lt increases. This is the case for all the models with 
f2/f2 cr it > 0.3 for Z > 0.00001. For Pop III stars, the increase of f/f C rit 
only occurs for the model with 0/J7 cr i t = 0.7. This is a consequence of the 
Gratton-Opick effect making the meridional circulation velocity smaller at 
lower Z. 



12 



3. At the end of the Main-Sequence phase, the increase of v /fcrit accelerates 
when the star contracts after the blue hook in the HR diagram. 

At solar metallicity, the lower initial value of S7/f2 cr it for a 20 M star to reach 
the critical limit during the Main- Sequence phase is about 0.8 (t>i n i =~ 400 km 
s" 1 ). At Z = 0.002 the limiting value is around 0.7. At Z=0.00001, the limiting 
value is below 0.7, while for Pop III stars, the limiting value is about 0.99 ! More 
precise values will be given in Ekstrom et al. (in preparation), however based on 
these results we can already draw the following consequences: only stars with 
a sufficiently high initial velocit y will reach the critic al limit. This is well in 
line with the results obtained bv lMartavan et aT. (2005) showing that the initial 



velocities of the Be stars is significantly higher than the initial velocities of the 
normal B stars. For 20 M Q stars, the range of initial velocities for which the 
star reaches the critical limit during the MS phase is larger at lower Z when Z 
varies between 0.02 and 0.00001. For pop III 20 M Q stellar models, this range 
is much smaller than at higher metallicities. For higher initial mass stars, the 
luminosity may be sufficiently close from t he Eddington limit that the fir-limit 
can be reached ([Maeder &: Mevnetll20 00bb 

The complete coverage of the ranges of initial masses, metallicities and 
velocities will allow to make predictions for the frequency of Be stars resulting 
from single star evolution at various metallicities. 
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